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ABSTRACT 


Mesoscale  Convective  Complex  vs.  Non-Mesoscale 
Convective  Complex  Thunderstorms: 

A  Comparison  of  Selected  Meteorological  Variables.  (August  1986) 

Michael  Eugene  Hoofard,  A.S.,  San  Antonio  College 
B.S.,  Texas  A4M  University 

Chairman  of  Advisory  Committee:  Prof.  Walter  K.  Henry 

/^r  A  comparative  investigation  of  mesoscale  convective  complex 
(MCC)  and  non-mesoscale  convective  complex  (non-MCC)  prestorm  environ¬ 
ments  is  conducted.  Eleven  atmospheric  variables  normally  associated 
with  thunderstorm  formation  are  either  observed  or  calculated  for  a 
total  of  nine  MCC  and  nine  non-MCC  storms.  These  variables  Include: 

850  mb  mixing  ratio,  850  mb  advection  of  water  vapor,  850  mb  flux 
divergence  of  water  vapor,  surface  to  500  mb  average  relative  humidity, 
precipitable  water,  the  Totals  Index,  the  Total  Energy  Index,  1000-700 
mb  thickness  advection  by  the  850  mb  wind,  500-300  mb  thickness  advec¬ 
tion  by  the  400  mb  wind,  vortlclty  advection  at  500  mb,  and  vertical 
velocity  at  700  mb.  Mean  values  are  calculated  for  the  eleven  varia¬ 
bles  according  to  storm  type.  Then,  the  mean  values  from  the  MCC  cases 
are  compared  statistically  with  the  corresponding  non-MCC  values. 

Results  show  that  a  significant  difference  exists  for  the  following 
mean  values:  850  mb  mixing  rat^o,  °^0  mb  advection  of  water  vapor, 
precipitable  water.  Total  Energy  Index,  and  1000-700  mb  thickness 
advection.  The  low-level  water  vapor  advection  and  thickness  advection  — v7 
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^Variables  are  combined  to  form  a  low-level  energy  rate  of  mange  term. 
This  energy  rate  of  change  term  Is  found  to  provide  an  even  better 
distinction  between  MCC  and  non-MCC  storm  environments. 
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1.  INTRODUCTION 

Prior  to  19/ j,  meteorologists  had  recognized  three  basic  types 
of  thunderstorm  systems.  The  classification  scheme  was  based  upon  the 
storm's  physical  appearance  and  organization  and  depended  largely  upon 
the  availability  of  radar  observations.  The  three  basic  types  are  the 
individual  air  mass  thunderstorm,  the  thunderstorm  cluster,  and  the 
squall  line.  Of  these  three,  the  squall  line  displayed  the  greatest 
organization  and  generally  covered  the  largest  area.  With  the  increas¬ 
ing  use  of  meteorological  satellites,  however,  a  fourth  type  has  been 
observed.  Referred  to  as  a  mesoscale  convective  complex  (MCC),  this 
storm  type  is  similar  to  the  squall  line  in  that  it  covers  large  areas 
and  appears  highly  organized.  However,  its  nearly  circular  shape  as 
depicted  in  satellite  imagery  readily  distinguishes  it  from  the 
elongated  squall  line. 

Primarily  nocturnal  in  nature,  MCCs  have  been  observed  over  the 
United  States  during  the  months  of  March  through  September.  Climato¬ 
logical  studies  for  the  United  States  (Maddox,  1981;  Maddox  et  al . , 
1982;  Rogers  and  Howard,  1983)  indicate  an  average  occurrence  of  33 
storms  per  season.  Of  the  131  MCCs  documented,  46  percent  produced 
tornadoes,  58  percent  produced  hail,  and  65  percent  produced  damaging 
winds.  In  addition,  flooding  occurred  with  36  percent  of  the  storms. 

Because  of  the  MCCs  destructive  potential,  interest  in  forecast¬ 
ing  these  intense  convective  storms  has  developed  quickly.  Unfortu¬ 
nately,  numerical  models  presently  in  operation  have  not  provided 

Citations  follow  the  style  of  the  Monthly  Weather  Review. 
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adequate  guidance  (Maddox,  1981).  Their  Inadequacy  should  not  be 
surprising,  however.  The  problem  plaguing  any  thunderstorm  forecast 
is  the  dependence  upon  synoptic  scale  data  to  predict  a  mesoscale 
event.  To  complicate  matters,  forecasters  attempt  to  predict  not  only 
thunderstorm  occurrence  but  the  level  of  severity  as  well.  Now,  with 
the  goal  of  forecasting  MCCs,  meteorologists  face  an  even  more  demand¬ 
ing  problem.  Can  the  atmospheric  environment  leading  to  MCC  formation 
be  distinguished  from  that  leading  to  other  types  of  convective  storms, 
severe  or  otherwise?  While  experiments  In  forecasting  MCCs  have  been 
undertaken,  no  one  has  addressed  this  question  directly.  This  Investi¬ 
gation  Is  a  beginning  venture  to  answer  this  question. 
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2.  LITERATURE  REVIEW  OF  THE  MCC 

MCCs  were  first  detected  and  recognized  as  unique  storm  systems 
during  the  mid  1970s  through  the  use  of  geosynchronus  satellite 
imagery.  Figure  1  depicts  a  typical  enhanced  infrared  (IR)  image  of 
the  MCC  while  Fig.  2  shows  the  corresponding  radar  echo.  The  formal 
definition  of  the  MCC  was  developed  solely  from  enhanced  IR  observa¬ 
tions  and  is  reproduced  in  Table  1. 

Maddox  (1981),  in  his  dissertation,  completed  the  first  compre¬ 
hensive  investigation  of  these  storms.  His  work  included  a  brief 
climatology  of  MCCs,  an  objective  analysis  of  the  meteorological 
conditions  spanning  their  life-cycle,  the  development  of  a  physical 
model  of  the  MCC,  and  an  introductory  investigation  of  the  storm's 
moisture  and  kinetic  energy  budgets.  He  found  the  following  features 
to  precede  persistently  the  genesis  of  these  convective  systems:  1) 
an  approaching  weak,  mid-tropospheric  short-wave  trough,  2)  strong 
low-level  warm  a dvection,  3)  a  conditionally  unstable  atmosphere,  4)  a 
significant  east-west  moisture  gradient,  and  5)  the  presence  of  a 
frontal  zone.  Additional  studies  by  Maddox  £t  aK  (1981)  dealt  with 
the  effects  of  the  MCC  on  its  environment.  They  concluded  that  the 
mesoscale,  convectively  driven  circulation  significantly  alters  upper- 
tropospheric  environmental  conditions,  particularly  temperature, 
height,  and  wind  fields. 

MCC  precipitation  studies  by  Fritsch  j?t£K  (1981)  showed  the 
beneficial  aspect  of  these  storm  systems.  MCCs  were  found  to  produce 
rainfall  over  large  areas  and  to  account  for  a  significant  portion  of 
the  rainfall  during  the  growing  season  over  much  of  the  corn  and  wheat 


Fig.  1.  Enhanced  infrared  satellite  image  for  1130  GMT  13  May 
1981.  MCC  identified  by  large,  black  region  centered  over  Missouri 
(after  Maddox  et  al . ,  1982). 


Table  1.  Definition  of  mesoscale  convective  complex  (after  Maddox, 
1981). 


Physical  characteristics 


Size:  A:  Cloud  shield  with  IR  temperature  ^  -32°C  must  have 

an  area  >  100,000  km2 

B:  Interior  cold  cloud  region  with  temperature  ^  -52°C 
must  have  an  area  >  50,000  km2 

Initiate:  Size  definitions  A  and  B  are  first  satisfied 

Duration:  Size  definitions  A  and  B  must  be  met  for  a  period  of  >6  h 

Maximum  Contiguous  cold  cloud  shield  (IR  temperature  <  -32°C) 
extent:  reaches  maximum  size 

Shape:  Eccentricity  (minor  axis/major  axis)  ^  0.7  at  time  of 

maximum  extent 

Terminate:  Size  definitions  A  and  B  no  longer  satisfied 


belts  of  the  United  States. 

Wetzel  ^t  aK  (1982)  investigated  the  dynamic  structure  of 
several  MCCs.  They  then  compared  the  MCCs'  vertical  profiles  of 
divergence  to  those  from  tropical  cloud  clusters  and  a  typical  mid¬ 
latitude  cyclone.  Finally,  they  compared  the  MCCs'  and  cyclone's 
vertical  profiles  of  average  meridional  sensible  heat  transport. 
Results  suggest  that  the  MCC  develops  in  a  highly  barotropic  environ¬ 
ment  and,  thus,  is  driven  basically  by  convective  instability. 

Cotton  £t  £l_.  (1983)  and  Wetzel  et  a±.  (1983),  using  data 
collected  during  the  South  Park  Area  Cumulus  Experiment,  investigated 
the  evolution  of  MCCs  from  mountain  generated  convective  cells  over 
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Colorado.  They  concluded  that  a  mid-latitude  MCC  is  similar  to  its 


& 


tropical  counterpart  in  that  their  dynamics  are  essentially  controlled 
by  buoyant  instability.  Thus,  they  suggested  that  the  weak  baroclinic 
zones  often  located  near  developing  MCCs  "act  primarily  to  trigger  and 
direct  the  release  of  buoyant  instability."  Cotton  et  al_.  further 
hypothesized  that  as  barocl inicity  increases,  convection  will  tend  to 
favor  the  typical  squall  line  structure  rather  than  the  elliptical 


V. 

['.< 


shape  characteristic  of  MCCs. 


During  the  summers  of  1982  and  1983,  experiments  in  forecasting 
the  MCC  were  conducted  jointly  by  personnel  at  the  Environmental 
Research  Laboratory  in  Boulder,  Colorado,  and  at  the  Satellite  Field 
Service  Station  in  Kansas  City,  Missouri.  A  summary  of  the  1982 
experiment  was  presented  at  the  American  Meteorological  Society's  13th 
Conference  on  Severe  Local  Storms  ( Maddox  ^t  al_. ,  1983).  During  this 
first  summer,  a  total  of  117  forecasts  were  made  with  the  following 
results:  20  forecasts  for  MCC  occurrence  verified,  58  forecasts  for 
no  occurrence  verified,  31  forecasts  for  MCC  occurrence  did  not  verify, 
and  8  MCC  events  were  not  forecast.  The  false  alarm  ratio,  the  number 
of  times  the  event  was  forecast  but  did  not  occur  divided  by  the  total 
number  of  times  the  event  was  forecast,  was  0.61.  Distinguishing 
differences  between  the  synoptic  conditions  leading  to  MCCs  and  those 
resulting  in  non-MCC  thunderstorms  proved  to  be  a  difficult  task. 
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3.  OBJECTIVE 

The  effects  of  the  MCC  on  man  and  hts  environment  have  been  both 
harmful  and  beneficial.  Consequently,  the  National  Oceanic  and 
Atmospheric  Administration  has  given  high  priority  to  the  development 
of  operational  MCC  forecast  procedures  and  techniques.  The  primary 
objective  of  this  research  Is  to  determine  if  MCC  and  non-MCC  prestorm 
environments  can  be  distinguished  from  one  another  using  routine 
synoptic  data.  To  accomplish  this  goal,  a  comparative  investigation 
of  MCC  to  non-MCC  convective  storms  was  performed. 
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4.  INVESTIGATIVE  PROCEDURE 

In  conducting  this  comparative  investigation,  selected  meteorolo¬ 
gical  variables  were  either  observed  or  calculated  for  a  total  of  nine  i 

i 

MCC  and  nine  non-MCC  storms.  Mean  values  were  then  determined  for  each  j 

i 

variable.  The  mean  values  derived  from  the  MCC  cases  were  compared  , 

i 

statistically  with  the  corresponding  non-MCC  values  to  determine  if  a 

significant  difference  existed.  The  following  sections  outline  this 

procedure  in  greater  detail.  j 

I 

a.  Selection  of  Meteorological  Variables  ! 

Thunderstorm  formation  depends  upon  the  availability  of 
sufficient  water  vapor,  unstable  air,  and  a  triggering  mechanism. 

These  three  requirements  guided  the  selection  of  the  meteorological 
variables.  A  summary  of  the  variables  selected  is  provided  in  Table  2. 

Table  2.  Meteorological  variables  selected  for  evaluation  of  MCC  and 
non-MCC  storms. 


1.  850  mb  mixing  ratio 

2.  850  mb  advection  of  water  vapor 

3.  850  mb  flux  divergence  of  water  vapor 

4.  Surface-to-500  mb  average  relative  humidity 

5.  Preci pi  table  water 

6.  Totals  Index 

7.  Total  Energy  Index 

8.  1000-700  mb  thickness  advection  by  850  mb  wind 

9.  500-300  mb  thickness  advection  by  400  mb  wind 

10.  Vorticity  advection  at  500  mb 

11.  Vertical  velocity  at  700  mb 


(1 )  Water  vapor 

By  mass,  approximately  one  half  of  the  atmosphere's  water  vapor 
is  located  below  800  mb.  Palmen  and  Newton  (1969)  have  estimated  that 
for  a  large  intense  thunderstorm,  about  90  percent  of  the  inflow  water 
vapor  enters  near  the  base  of  the  updraft.  Sienkiewicz  (1981),  in  her 
moisture  budget  study  of  convective  activity,  "found  horizontal 
moisture  convergence  to  be  the  dominant  term  within  the  900-750  mb 
layer."  Maddox  (1981)  observed  similar  results  in  his  investigation 
of  ten  MCCs .  Consequently,  low-level  water  vapor  measurements  and 
processes  have  been  stressed  in  this  study.  The  850  mb  mixing  ratio 
observations  were  used  to  obtain  an  instantaneous  measure  of  available 
water  vapor;  850  mb  water  vapor  advection  and  flux  divergence  were 
calculated  to  determine  changes  occurring  within  the  moisture  field. 

A  question  raised  during  this  research  concerned  the  adequacy  of 
the  850  mb  mixing  ratio  in  representing  low-level  water  vapor  content. 
More  specifically,  how  often  might  a  moist  surface  layer  go  unrecog¬ 
nized  due  to  the  capping  effect  of  an  inversion  below  the  850  mb 
surface?  The  mean  mixing  ratio  of  the  lowest  100  mb  is  a  more  repre¬ 
sentative  value  of  low-level  water  vapor  content  but  it  is  not  readily 
available.  To  determine  if  the  850  mb  data  adequately  depicted  the 
presence  of  low-level  humidity,  a  comparison  was  made.  Dewpoint 
depressions  at  the  850  mb  level  were  compared  against  the  corresponding 
surface  values  for  each  storm  case  selected.  For  the  MCC  data,  16 
percent  of  the  observations  had  850  mb  dewpoint  depressions  more  than 
3°C  larger  than  the  corresponding  surface  values.  For  the  non-MCC 
data,  a  value  of  13  percent  was  found.  In  both  cases.  Gulf  coastal 


stations  accounted  for  approximately  45  percent  of  the  dewpoint 
depression  variations.  None  of  the  storm  cases  used  in  this  research, 
however,  occurred  within  550  km  of  the  coast.  Based  upon  these 
findings,  the  850  mb  mixing  ratio  was  accepted  as  a  sufficiently 
accurate  Indicator  of  low-level  humidity. 

Surface-to-500  mb  average  relative  humidity  and  preci pi  table 
water  were  selected  as  additional  Indicators  of  available  water 
vapor.  These  two,  however.  Include  the  added  effect  of  middle  and 
upper-level  water  vapor.  Maddox  (1983)  found  MCCs  to  form  near 
preci pi  table  water  maxima. 

(2)  Stability 

Both  conditional  and  convective  stabilities  play  Important  roles 
in  the  promotion  of  convection.  If  the  environmental  lapse  rate  lies 
between  the  dry  and  moist  adiabatic  lapse  rates,  the  atmosphere  is 
said  to  be  conditionally  unstable.  With  sufficient  lifting,  a  parcel 
of  air  becomes  saturated  and  eventually  buoyant  with  respect  to  its 
surroundings.  Often,  however,  soundings  which  are  conditionally 
unstable  in  the  lower  levels  may  exhibit  capping  Inversions  at  higher 
levels.  Unless  the  Inversion  is  removed,  convection  will  be 
Inhibited.  If  the  equivalent  potential  temperature  within  the 
Inversion  layer  decreases  with  height,  upward  displacement  of  the 
layer  will  result  in  destabil ization  and,  thus,  destroy  the  inversion. 
Layers  of  the  atmosphere  In  which  the  equivalent  potential  temperature 
decreases  with  height  are  said  to  be  convectively  unstable. 

Various  Indices  have  been  developed  to  depict  the  degree  of 
atmospheric  Instability.  In  general,  these  Indices  yield  similar 
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results,  especially  for  a  warm,  moist  air  mass  (Miller,  1967).  Two 
indices  selected  for  this  study  are  the  Totals  Index  and  Total  Energy 
Index.  A  study  of  the  1964  and  1965  tornado  events  revealed  that  92 
percent  of  the  storms  were  characterized  by  a  Totals  Index  of  50  or 
greater  (Miller,  1967).  The  Total  Energy  Index,  which  is  based  upon 
the  vertical  variation  in  static  energy,  has  also  proved  useful  in 
identifying  severe  storm  threat  areas.  Darkow  (1968)  found  total 
energy  values  to  be  significantly  higher  for  soundings  near  tornadic 
storms.  In  addition,  he  maintains  that  the  index  holds  a  slight 
advantage  over  the  Showalter  and  Lifted  indices  through  its  ability  to 
account  for  contributions  of  cold  mid-tropospheric  air.  Intrustions 
of  cold,  upper-level  air  appear  to  play  an  important  role  in  the  total 
energy  release  of  severe  convective  storms.  Tables  3  and  4  show  the 
guidance  provided  for  applying  the  Totals  and  Total  Energy  indices. 

(3)  Triggering  actions 

The  presence  of  abundant  water  vapor  and  potentially  unstable 
layers  does  not  ensure  the  development  of  thunderstorms.  Unless  some 
mechanical  process  induces  upward  vertical  motion,  the  potential 
instability  will  not  be  released.  Triggering  actions  evaluated  in  this 
study  include  low-level  warm  advection,  upper-level  cold  advection,  and 
positive  vorticity  advection. 

The  quasi -geostrophic  omega  equation  is  useful  for  diagnosing 
vertical  motion  fields  within  the  large  scale  baroclinic  environment. 
Qualitatively,  air  rises  in  regions  where  either  differential  vorticity 
advection  or  low-level  warm  advection  is  occurring  (Holton,  1972). 
Quantitative  assessment  is  difficult,  however,  since  the  differential 


Table  3.  Suggested  guidance  for  using  the  Totals  Index  (after  Miller, 
1967). 

i 

I 


Totals  Index  Value 

Interpretation 

44 

isolated  thunderstorms 

46 

scattered  thunderstorms;  some  heavy 

48 

scattered  thunderstorms;  some  heavy  and 
isolated  severe 

50 

scattered  heavy  thunderstorms,  some  severe 

52 

scattered  to  numerous  heavy  thunderstorms; 
some  severe  with  tornadoes  possible 

56 

numerous  heavy  thunderstorms;  scattered 
severe  thunderstorms  with  tornadoes 

Table  4.  Suggested  guidance  for  using  the  Total  Energy  Index  (after 
Darkow,  1968). 


Total  Energy  Index  Value  Interpretation 

0.0  to  -1.0 
-1.0  to  -2.0 


less  than  -2.0 


thunderstorms  possible;  not  severe 
thunderstorms;  isolated  severe 
severe  thunderstorms;  tornadoes  likely 
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vorticity  advection  term  and  temperature  advection  term  often  imply 
vertical  motions  of  opposite  sign.  Studies  of  the  MCC  prestorm 
environment  (Maddox,  1983)  generally  indicate  the  presence  of  strong 
warm  advection  below  700  mb  but  quite  weak  PVA  aloft.  These  findings 
suggest  that  the  vertical  motion  preceding  MCC  development  is  induced 
primarily  by  warm  advection. 

The  validity  of  applying  the  omega  equation  to  the  MCC  environ¬ 
ment  is  questionable,  however.  Foremost,  findings  from  earlier 
investigations  tend  to  classify  the  MCC  environment  as  barotropic  in 
spite  of  the  storm's  usual  proximity  to  a  surface  front  and  advancing 
short-wave  trough.  But  even  if  the  environment  can  be  classified  as 
baroclinic,  the  proximity  to  a  surface  front  raises  scaling  questions. 
Derived  for  the  diagnosis  of  mid-latitude  synoptic  scale  motions,  the 
omega  equation  does  not  include  the  vertical  advection  and  twisting 
terms  of  vorticity.  These  terms  are  often  important  in  the  vicinity 
of  atmospheric  fronts  (Holton,  1972).  Thus,  because  the  MCC  prestorm 
environment  often  includes  a  surface  front,  the  omega  equation  may  not 
be  applicable. 

The  quantitative  relationship  between  low-level  warm  advection 
and  the  vertical  motion  field  associated  with  developing  MCCs  is  still 
unknown.  The  destabilizing  effect  of  warm  advection  and  its  relation¬ 
ship  to  severe  convective  activity  has  been  observed  for  many  years, 
however.  As  early  as  1952,  MacDonald  related  lower- tropospheric  warm 
advection  to  the  formation  and  maintenance  of  squall  lines.  Recently, 
Hales  (1982)  demonstrated  the  predominance  of  strong  warm  advection  at 
both  850  and  700  mb  during  to.  nado  events.  The  close  relationship 
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which  has  been  observed  between  warm  advection  and  severe,  organized 
convection  warrants  its  selection  as  a  possible  triggering  action. 

Since  cold  advection  aloft  also  results  in  atmospheric 
destabilization,  it  too  will  be  evaluated  as  a  triggering  action. 
Crawford  (1950)  noted  the  common  occurrence  of  cold  advection  aloft 
during  the  formation  of  pre-frontal  thunderstorms.  He  demonstrated 
cases  of  differential  advection  in  which  the  cold,  upper-level  trough 
moved  eastward  faster  than  the  low-level  warm  sector.  Fulks  (1951) 
made  similar  observations  in  his  study  of  the  instability  line.  Rhea 
(1966)  noted  that  thunderstorms  often  formed  along  dry  lines  when 
supported  by  cold  advection  aloft.  Since  cold  pockets  aloft  are 
usually  associated  with  troughs  and  minor  short-waves,  it  may  be 
technically  more  appropriate  to  classify  cold  advection  as  an  Indicator 
rather  than  an  actual  trigger  (Palmen  and  Newton,  1969).  Thus  a  dynam¬ 
ical  triggering  process  needs  to  be  considered  as  well. 

Basic  dynamic  theory  depicts  rising  air  downstream  of  upper-level 
troughs  and  descending  air  downstream  of  upper-level  ridges.  The 
presence  of  these  vertical  motion  fields  is  often  explained  in  terms 
of  vorticity  advection.  On  the  synoptic  scale,  PVA  downstream  of  the 
trough  generally  implies  the  presence  of  convergence  below  the  level 
of  nondivergence  and  divergence  above.  Thus,  the  upward  vertical 
motion  field  is  established.  Miller  (1967)  ranked  vorticity  advection 
as  the  most  important  variable  related  to  severe  weather  out-breaks. 
Although  more  recent  studies  (Hales,  1979;  Maddox  and  Doswell,  1982) 
downplay  the  importance  of  PVA  at  500  mb,  it  was  selected  as  a  trigger¬ 
ing  action  to  be  evaluated. 
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(4)  Vertical  motion 

Sustained  upward  motion  Is  critical  to  the  formation  and  persist¬ 
ence  of  thunderstorm  activity.  In  general,  the  stronger  the  upward 
velocity,  the  more  favorable  an  area  will  be  for  thunderstorm  forma¬ 
tion.  Vertical  velocity  was  the  final  meteorological  variable  examined 
for  its  value  in  distinguishing  between  MCC  and  non-MCC  prestorm 
environments. 

b.  Data  Sources 

Values  for  the  selected  variables  were  determined  using  either 
National  Meteorological  Center  (NMC)  facsimile  charts  or  the  mandatory 
level  sounding  data.  Precl pi  table  water  and  surface-to-500  mb  average 
relative  humidity  data  came  directly  from  the  NMC  charts.  Vertical 
velocity  was  estimated  using  the  LFM  12-24  h  forecast  charts. 

Vorticity  advection  was  estimated  using  the  500  mb  Helght/Vortlcity 
Analysis  chart.  The  remainder  of  the  variables  were  either  obtained 
directly  or  computed  from  the  sounding  data.  Storm  locations  were 
identified  using  radar  facsimile  charts.  For  informational  purposes, 
NMC  surface  charts  were  reviewed  to  identify  frontal  locations. 

c.  Case  Selection  Criteria 

Because  sounding  data  were  available  only  for  those  stations 
identified  in  Fig.  3,  storms  occurring  outside  this  domain  or  along 
its  borders  were  not  considered.  MCC  and  non-MCC  selection  criteria 
differed  somewhat;  MCC  selection  is  discussed  first.  Only  those  MCCs 
identified  in  the  1981  and  1982  annual  summaries  (Maddox  et  al.,  1982; 


Rogers  and  Howard,  1983)  were  considered.  Thus,  the  storms  selected 
as  MCCs  had  already  been  verified  as  meeting  the  formal  criteria  used 
to  define  the  MCC.  These  summaries  listed  both  the  time  in  which  the 
first  storms  occurred  and  the  time  in  which  the  system  was  classified 
as  an  MCC.  Only  those  systems  with  first  storms  occurring  after  0000 
GMT  but  before  1200  GMT  were  selected.  As  noted  in  the  literature 
review,  the  MCC  has  been  found  to  alter  the  surrounding  upper- 
tropospheric  environmental  conditions.  Since  this  investigation 
emphasizes  the  pre-storm  environment,  an  attempt  was  made  to  exclude 
or  at  least  minimize  possible  storm  impacts.  The  annual  summaries 
indicate  that  approximately  75  percent  of  the  MCCs  develop  between 
0000  and  1200  GMT;  over  50  percent  develop  between  0000  and  0600  GMT. 
Using  0000  GMT  sounding  data  and  selecting  storms  occurring  after  this 
time  maximized  the  number  of  cases  available  in  the  late  prestorm 
stage.  To  limit  the  amount  of  elapsed  time  between  the  0000  GMT 
sounding  and  the  time  of  MCC  formation,  no  MCCs  forming  later  than 
1200  GMT  were  considered;  1200  GMT  sounding  data  were  avoided  due  to 
the  greater  possibility  of  low-level  inversions  and  unrepresentative 
humidity  values  at  850  mb.  The  1200  GMT  data  were  used  in  one  case. 
For  the  18  July  1982  MCC,  the  preceding  1200  GMT  surf ace- to-500  mb 
relative  humidity  and  precipitable  water  charts  were  substituted  in 
place  of  the  missing  0000  GMT  data.  Comparison  of  the  substitute 
chart  with  that  available  24  hours  later  showed  little  change  of 
patterns.  Thus,  the  morning  sounding  prior  to  the  storm  was  deemed  a 
representative  replacement.  Ten  cases  were  identified  using  these 
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criteria;  however,  one  had  to  be  omitted  due  to  missing  data.  Table  5 
lists  the  nine  MCC  events  selected. 


Table  5.  Times  and  dates  of  storms  investigated. 


MCC  Cases 

Non-MCC  Cases 

1. 

0735  GMT  13  May  1981 

0735  GMT  16  May  1979 

2. 

0835  GMT  7  June  1982 

0635  GMT  30  May  1979 

3. 

1035  GMT  8  June  1982 

0635  GMT  8  June  1978 

4. 

0635  GMT  10  June  1981 

0735  GMT  24  June  1979 

5. 

1035  GMT  22  June  1981 

0635  GMT  11  July  1978 

6. 

0635  GMT  18  July  1982 

0735  GMT  16  July  1979 

7. 

0935  GMT  23  July  1981 

0835  GMT  31  July  1978 

8. 

0835  GMT  15  August  1982 

0435  GMT  9  August  1978 

9. 

0535  GMT  2  September  1982 

0835  GMT  2  September  1978 

Selection  of  non-MCC  cases  was  more  complicated.  To  ensure  that 
the  non-MCC  storms  were  independent  of  the  MCC  cases  selected,  1981 
and  1982  storm  data  were  not  considered.  However,  the  problem  of 
separating  non-MCC  storm  environments  still  remained.  To  resolve  this 
problem,  MCC  annual  summaries  (Maddox,  1981)  for  1978  and  1979  were 
reviewed.  For  this  two-year  period,  an  attempt  was  made  to  select 
non-MCC  storms  occurring  at  least  24  h  before  or  after  an  MCC  event. 

All  except  one  of  the  cases  eventually  chosen  met  this  requirement. 

The  exception  occurred  approximately  21  h  after  an  MCC  but  was  included 
to  increase  the  number  of  non-MCC  cases.  Since  0000  GMT  data  were  used 
for  the  non-MCC  cases  as  well,  selection  of  storms  forming  between 
0000-1200  GMT  was  preferred.  However,  in  order  to  obtain  a  sufficient 
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number  of  non-MCC  cases,  storms  which  formed  prior  to  0000  GMT  but 
showed  development  and  Intensification  during  the  preferred  time 
period  were  accepted.  Table  5  also  lists  the  nine  non-MCC  cases 
selected. 

d.  Evaluation  of  Computed  Variables 

For  each  case,  sounding  data  from  the  rawlnsonde  station  network 
were  entered  Into  a  computer.  Corresponding  850  and  500  mb  mixing 
ratio,  surface-to-500  mb  relative  humidity,  and  precipltable  water 
data  had  to  be  entered  as  well.  Thickness  values  were  determined  from 
the  raw  data  for  both  the  1000-700  and  500-300  mb  layers.  In  addition, 
the  850  and  400  mb  wind  speeds  were  converted  from  kt  to  km  h~^  and 
separated  Into  u  and  v  components. 

The  Totals  Index  (TI)  was  calculated  for  each  station  as  shown 
by 


TI  =  (TD85  -  T50)  +  (T85  -  T50)  (1) 

»  TD85  +  T85  -  2  T50 

where  TDQ5>  Tg5,  and  TgQ  are  the  850  mb  dewpoint,  850  mb 
temperature,  and  500  mb  temperature  In  °C. 

The  Total  Energy  Index  was  calculated  using  the  two-step  procedure 
shown  In  (2)  and  (3).  First,  the  total  static  energy  (TE)  per  unit 
mass  was  determined  for  both  the  850  and  500  mb  levels  by 


TE  =  . 24( T  ♦  2.5w  ♦  9.8z). 


(2) 
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For  the  specified  pressure  level,  T  is  the  temperature  in  °K,  w  Is 
the  mixing  ratio  in  g  kg~\  and  z  is  the  height  in  km.  Units  for  TE 
are  cal  g"1.  See  Oarkow  (1968)  for  details  concerning  the  develop¬ 
ment  of  this  equation.  By  subtracting  the  850  mb  static  energy  total 
from  the  500  mb  value,  the  Total  Energy  Index  is  obtained. 

TEI  =  TE50  -  TE85  (3 

At  this  point,  Barnes'  objective  analysis  (Barnes,  1964)  was 
applied  to  interpolate  the  following  station  data  to  a  grid:  1)  850 
mb  mixing  ratio,  2)  surface-to-500  mb  relative  humidity,  3)  precipi- 
table  water,  4)  Totals  Index,  5)  Total  Energy  Index,  6)  u  and  v 
components  of  the  850  and  400  mb  wind,  and  7)  thickness  values  for  the 
100-700  and  500-300  mb  layers.  The  grid  was  designed  for  use  with  the 
NMC  radar  facsimile  chart.  This  feature  allowed  direct  comparison  of 
analyzed  fields  to  the  associated  storm  echoes.  The  grid  is  illus¬ 
trated  in  Fig.  4;  grid  spacing  is  250  km  at  40°  north  latitude. 

Upon  completion  of  the  interpolation  scheme,  grid  values  had  been 
determined  for  five  of  the  eleven  selected  variables.  The  gridding 
technique  was  not  used  in  estimating  vorticity  advection  or  vertical 
motion.  However,  grid  values  were  required  for  these  remaining 
variables:  1)  low-level  warm  advection,  2)  upper-level  cold 
advection,  3)  850  mb  advection  of  water  vapor,  and  4)  850  mb  flux 
divergence  of  water  vapor.  The  computations  of  these  variables  are 
discussed  below. 

Temperature  advection  in  the  desired  layers  was  determined  in 


terms  of  thickness  advection.  Low-level  thickness  advection  was 
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estimated  using  the  1000-700  mb  thickness  and  the  850  mb  wind.  Upper- 
level  thickness  advection  was  estimated  using  the  500-300  mb  thickness 
and  the  400  mb  wind.  Thickness  advection  values  were  calculated  for 
the  center  of  each  grid  box.  Advection  by  the  u  and  v  wind  components 
was  calculated  separately  and  then  added  to  obtain  the  total  advection. 
Equations  (4),  (5),  and  (6)  In  combination  with  Fig.  5  illustrate  the 
computational  method  employed. 


Fig.  5.  Sample  grid  box  with  location  of  variables  used  for 
thickness  advection  calculation. 
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Thickness  advection  by  u: 

-u  ( a  a  0)  _  f-(u]  +  U2  +  U3  +  U4 )  1  T(a$2  +  A04 )  -  (a0i  ♦  a03) 

ax  1  4  J  L  2ax 

Thickness  advection  by  v: 

-v  (j>A0)  =  [-(vi  +  vg  ♦  v3  +  V4 )1  [(a#i  ♦  4^2)  -  (a03  ♦  a04) 

ay  l  4  J  [  2Ay 

Total  advection: 

(A)  (B) 

-V-v(a 0)  =  - [u( a  A0)  +  v( a A0) 1  (6) 

I  ax  ay  J 

where  terms  A  and  B  represent  the  values  of  expressions  (4)  and  (5), 
respectively.  For  these  and  the  following  calculations,  units  for  x 
and  y,  u  and  v,  and  0  are  km,  km  h'1 ,  and  gpm.  Resulting  units  for 
thickness  advection  are  gpm  h"1. 

Computation  of  the  advection  of  water  vapor  at  850  mb  was  per¬ 
formed  in  a  similar  manner.  The  substitution  of  mixing  ratio  (w)  for 
thickness  (a0)  was  the  only  change  required  before  using  (4),  (5),  and 
(6).  Units  are  g  kg-1  for  w  and  g  kg"1  h"1  for  the  advection  of  water 
vapor. 

Flux  divergence  of  water  vapor  at  850  mb  was  calculated  using  (7). 
Like  the  advection  values,  it  was  determined  for  the  center  of  each 
grid  box.  For  all  computations  involving  water  vapor,  the  mixing 
ratio  has  been  substituted  in  place  of  the  specific  humidity  (q). 
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w?*V  =  (wi+W2+W3+W4)  f{u2+U4)-(u]+U3)  ♦  (vi+V2)-(v3+V4) 

4  l  2ax  2ay 

Resulting  units  for  the  flux  divergence  of  water  vapor  are  g  kg-1 
h  .  The  total  flux  divergence,  v*(wV),  was  not  included  among  the 
variables  to  be  investigated.  Maddox's  (1981)  investigation  of  ten  MCCs 
showed  little  numerical  difference  between  vapor  flux  divergence  and 
total  vapor  flux  divergence  in  the  MCC  prestorm  environment. 

e.  Evaluation  of  Gridded  Fields 

Using  the  processing  and  computational  procedures  discussed 
above,  grid  values  were  computer  produced  for  the  first  nine  variables 
listed  in  Table  2.  Thus,  nine  fields  were  obtained  for  each  MCC  and 
non-MCC  case.  Although  each  field  was  drawn  by  hand,  the  exact  loca¬ 
tion  of  the  associated  storm  was  not  identified  before  completing  the 
analysis.  This  technique  was  used  to  promote  analysis  objectivity. 

To  determine  the  variables  associated  with  each  storm,  each 
analyzed  chart  was  positioned  over  the  corresponding  radar  facsimile 
chart.  Radar  echoes  better  depicted  the  MCC's  region  of  strong, 
active  convection  and,  thus,  were  chosen  over  enhanced  IR  satellite 
imagery.  While  each  analyzed  chart  was  based  on  0000  GMT  data,  the 
corresponding  radar  chart  was  based  on  a  radar  observation  taken  some¬ 
time  between  0000-1200  GMT.  For  each  MCC  case,  the  actual  observation 
selected  was  the  one  taken  closest  to  the  time  of  MCC  formation.  As 
mentioned  earlier,  the  formation  times  were  available  within  the  MCC 
annual  summaries.  In  general,  the  radar  chart  selected  for  a  non-MCC 


case  was  one  depicting  storm  development  and  intensification.  A 
series  of  charts  for  the  MCC  and  non-MCC  conditions  are  shown  in 
Appendix  A.  Radar  echoes  not  a  part  of  the  MCC  or  non-MCC  storm 
system  are  not  shown. 

After  superimposing  each  set  of  isopleths  over  the  correspond¬ 
ing  storm  radar  echo,  representative  values  were  determined  using  the 
following  rules. 

1)  Only  values  occurring  within  the  echo  region  were  considered. 

2)  The  largest  850  mb  mixing  ratio,  surface-to-500  mb  relative 
humidity,  and  preci pi  table  water  values  were  recorded. 

3)  The  largest  positive  value  of  850  mb  water  vapor  advection 
was  recorded;  if  the  echo  region  contained  no  positive  values,  the 
smallest  advection  value  was  recorded. 

4)  The  largest  value  of  850  mb  water  vapor  flux  convergence  was 
recorded;  if  the  echo  region  contained  no  convergence  values,  the 
smallest  value  of  water  vapor  flux  divergence  was  recorded. 

5)  The  Totals  Index  and  Total  Energy  Index  values  representing 
the  greatest  instability  were  recorded. 

6)  The  largest  positive  value  of  low-level  thickness  advection 
was  recorded;  if  the  echo  region  contained  no  positive  values,  the 
smallest  advection  value  was  recorded. 

7)  The  largest  negative  value  of  upper-level  thickness  advection 
was  recorded;  if  the  echo  region  contained  no  negative  values,  the 
smallest  advection  value  was  recorded. 

Certain  limitations  and  biases  were  inherent  to  this  evaluation 
procedure.  Time  differences  between  the  analyzed  fields  and  the 
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corresponding  radar  echoes  limited  the  representativeness  of  the 
recorded  values.  This  limitation  was  not  considered  a  major  problem, 
however,  because  of  generally  light  mid-level  steering  winds  and  a 
relatively  large  grid  spacing  (250  km).  In  other  words,  storm  move¬ 
ment  between  0000  GMT  and  the  time  of  the  radar  depiction  was  normally 
less  than  one  grid  space.  The  procedure  for  selecting  largest  and 
smallest  values  produced  biased  results.  The  values  favorable  to 
thunderstorm  development  were  maximized  while  those  unfavorable  were 
minimized.  Since  this  technique  was  applied  to  both  MCC  and  non-MCC 
cases,  the  bias  was  not  considered  detrimental  to  the  investigation. 

The  flexibility  of  this  technique  helped  compensate  for  the  time 
limitation  discussed  above  and,  thus,  aligned  results  more  closely 
with  theory.  The  use  of  extreme  rather  than  average  values  was 
important  for  another  reason  as  well.  The  identification  of  a  meso- 
scale  disturbance  within  a  synoptic  network  is  more  likely  when  using 
data  extremes  instead  of  averages. 

f .  Estimation  of  Vorticity  Advection 

Vorticity  advection  estimates  were  obtained  for  each  case  using 
0000  GMT  NMC  500  mb  Height/Vorticity  Analysis  charts.  After  enlarging 
the  charts,  they  were  altered  to  include  30  m  contour  intervals.  The 
MCC  and  non-MCC  echoes  as  depicted  on  the  radar  charts  were  transferred 
to  the  corresponding  height/vorticity  charts.  E :ho  selection  times 
were  the  same  as  those  discussed  in  section  4e.  On  each  NMC  analysis, 
a  square  box  with  sides  of  500  km  (2  grid  spaces)  was  drawn  upstream 
from  the  radar  echo  as  indicated  by  the  contours.  The  storm  center 
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Fig.  6.  Illustration  of  technique  for  estimating  vorticity 
advection.  Solid  lines  are  500  mb  height  contours  in  gpm.  Dashed 
lines  are  isolines  of  vorticity  { X10~^s_l ) .  Point  A  represents 
the  storm  center.  Height  contour  and  vorticity  isoline  intersections 
within  the  box  are  highlighted  by  circles.  For  this  example,  the  boxed 
region  is  characterized  by  PVA  with  four  intersections. 


was  located  at  the  midpoint  of  the  box's  downstream  end.  This  down¬ 
stream  end  was  positioned  normal  to  the  contours.  Figure  6  illustrates 
such  a  box.  The  area  within  each  box  was  then  classified  according  to 
the  type  of  vorticity  advection.  Thus,  the  boxes  represented  areas  of 
either  positive,  negative,  or  neutral  vorticity  advection.  The 
intensity  of  the  PVA  or  NVA  was  evaluated  by  counting  the  number  of 
height  contour  and  vorticity  isoline  intersections.  A  larger  number 
of  intersections  implies  stronger  PVA  or  NVA.  This  procedure  provided 


a  quick  field  method  for  evaluating  the  type  and  intensity  of  vorticity 
advection  associated  with  each  storm. 

g.  Estimation  of  Vertical  Motion 

The  LFM  provides  12-48  h  vertical  velocity  forecasts  for  the  700 
mb  level.  Isolines  of  vertical  motion  are  drawn  in  intervals  of  2  ub 
s"1  beginning  with  the  +  0  line.  The  12  and  24  h  forecasts  valid  at 
0000  and  1200  GMT  were  used  to  infer  synoptic  scale  vertical  motions 
associated  with  each  storm.  A  single  value  for  vertical  velocity  was 
interpolated  using  both  charts  since  the  time  of  storm  occurrence 
varied  between  0000-1200  GMT.  Estimates  were  made  to  the  nearest 
ub  s-1. 

h.  Statistical  Evaluation 

Once  the  variables  had  been  obtained  for  each  storm,  mean  values 
were  calculated  for  the  MCC  and  non-MCC  cases.  Then,  the  MCC  mean 
values  were  compared  statistically  against  their  non-MCC  counterparts 
to  determine  if  a  significant  difference  existed.  Assuming  indepen¬ 
dence  and  equal  but  unknown  variances,  the  hypothesis  that  the  corre¬ 
sponding  MCC  and  non-MCC  mean  values  were  not  significantly  different 
was  tested  using  the  Student  t  test.  The  test  statistic  for  testing 
the  equality  of  means  is  shown  by 

t  =  (Xi  -  3?2)/[s2(l/n1  +  l/n2)]1/2  (8) 

where  and  X2  are  the  mean  values  for  each  test  group,  S^  is 
the  estimate  of  the  common  variance,  and  n^  and  n£  are  the  number 
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of  observations  from  each  group.  S  is  calculated  using 

=  [(nj  -  l)Si^  +  (ri2  -  l)S2^]/(nj  +  n2  -  2).  (9) 

Again,  nj  and  nz  are  the  number  of  observations  from  each  test 
group  while  and  S2^  are  the  sample  variances  for  each  test 
group. 

Once  the  test  statistic  for  each  pair  of  mean  values  was  com¬ 
puted,  it  was  checked  against  the  appropriate  critical  region  found  in 
a  cumulative  t  distribution  table.  If  the  test  statistic  was  located 
within  the  acceptance  region  for  a  specified  significance  level,  the 
hypothesis  of  equality  was  accepted  and  the  two  means  were  considered 
to  come  from  the  same  population.  If  the  test  statistic  was  located 
outside  of  this  region,  however,  the  hypothesis  was  rejected  and  the 
two  means  were  considered  to  come  from  different  populations.  When 
testing  the  equality  of  two  means,  the  acceptance  region  is  defined 
by  the  Interval 


±  t  ( l-a/2) (ni*n2-2)  (10) 

where  a  represents  the  significance  level.  The  significance  level  is 
the  probability  of  rejecting  a  hypothesis  that  is  actually  correct. 

For  this  investigation,  however,  a  particular  significance  level  was 
not  specified.  Instead,  the  significance  level  at  which  rejection 
first  occurred  was  recorded.  Only  levels  between  0  and  .50  were 
considered. 

This  statistical  analysis  was  based  on  a  total  of  nine  MCCs  and 
nine  non-MCCs.  The  small  sample  size  is  recognized  as  an  investigative 
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weakness  since  statistical  reliability  Is  lessened.  Currently,  a 
relaxation  In  case  selection  criteria  would  be  required  to  Increase  the 
MCC  sample  size.  As  more  MCC  annual  summaries  become  available,  the 
sample  size  problem  will  diminish. 


5. 


ERROR  DISCUSSION 


Fuelberg  (1974)  has  estimated  the  root  mean  square  (RMS)  errors 
associated  with  upper-air  measurements.  His  results  are  shown  in 
Table  6.  Based  upon  his  data,  errors  have  been  estimated  for  the 
observed  and  computed  variables  used  In  this  study.  Table  7  lists  the 
RMS  errors  determined  for  each  parameter.  Fuel  berg's  estimate  of 
humidity  error  was  applied  directly  to  the  first  three  variables. 

The  errors  for  the  stability  indices  were  estimated  using  the 
following  procedure.  Values  from  (1)  and  (3)  were  calculated  using 
typical  estimates  for  the  right-hand  side  quantities.  Thus,  a  typical 
value  was  obtained  for  the  two  indices.  The  right-hand  side  quantities 
were  then  altered  by  an  amount  equal  to  the  respective  RMS  errors  found 
in  Table  6.  Care  was  taken  at  this  step  to  ensure  that  the  inserted 
errors  would  be  additive  rather  than  compensating.  Using  these 
adjusted  quantities,  the  values  were  recalculated.  The  difference 
between  the  first  and  second  solution  for  each  index  represented  an 
estimate  of  the  RMS  error. 

RMS  errors  for  advectlon  and  flux  divergence  computations  were 
estimated  in  basically  the  same  manner.  Using  equations  (4)  through 
(7)  and  grid  point  values  from  data  fields  obtained  in  this  research, 
advection  and  divergence  were  computed.  After  this  calculation,  all 
of  the  quantities  at  a  single  grid  point  were  altered  by  an  amount 
equal  to  their  respective  RMS  errors.  Again,  care  was  taken  to  ensure 
the  inserted  errors  would  be  additive.  The  solution  was  then 
recalculated.  As  before,  the  difference  from  the  original  solution 
represented  the  error  estimate.  This  procedure  was  conducted 
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Table  6.  Estimated  root  mean  square  errors  associated  with  upper-air 
measurements  (after  Fuelberg,  1974). 


Parameter 


Pressure  Level  Approximate  RMS  Error 


Temperature 

Pressure 

Humidity 

Pressure  altitude 

i 

|  Wind  speed 


Wind  direction 


— 

0.5OC 

to  1 

400  mb 

1 

.3  mb 

to 

100 

mb 

1 

.1  mb 

to 

10 

mb 

0 

i.7  mb 

— 

10% 

500 

mb 

10 

'  gpm 

300 

mb 

20 

gpm 

50 

mb 

50 

gpm 

700 

mb 

0.5 

m  s~l 

for 

40° 

elevation 

2.5 

m  s_l 

for 

10° 

elevation 

500 

mb 

0.8 

m  s-1 

for 

40° 

elevation 

4.5 

m  s-1 

for 

10° 

elevation 

300 

mb 

1.1 

m  s_1 

for 

40° 

elevation 

7.8 

m  s~l 

for 

10° 

elevation 

700 

mb 

1.30 

for 

40° 

elevation 

9.50 

for 

10O 

elevation 

500 

mb 

1.80 

for 

40° 

elevation 

13.40 

for 

10O 

elevation 

300 

mb 

2.50 

for 

400 

elevation 

18.00 

for 

IOO 

elevation 

A'*.  •**.>**'■<•  i 
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Table  7.  Estimated  root  mean  square  errors  for  selected  variables. 


850  mb  mixing  ratio 

Surface-to-500  mb  average  relative  humidity 
Preci pi  table  water 
Totals  Index 
Total  Energy  Index 

850  mb  water  vapor  advection 

850  mb  flux  divergence  of  water  vapor 

1000-700  mb  thickness  advection 


500-300  mb  thickness  advection 


1.0  cal  g-1 
0.13  g  kg-*  h“l 
0.05  g  kg-1  h~l 
0.6  gpm  h“l 
3.0  gpm  h_1 


separately  for  all  four  grid  points  to  determine  the  maximum  error. 

The  entire  procedure  was  then  repeated  for  several  grid  boxes  to 
obtain  an  average  maximum  error. 

Table  7  represents  only  those  errors  caused  by  the  accuracy 
limitations  of  meteorological  instruments.  Errors  introduced  through 
the  use  of  the  gridding  scheme  were  not  considered.  Regardless,  the 
RMS  error  estimates  provided  should  still  be  useful  when  evaluating 
the  reliability  of  investigative  results.  Because  vorticity  advection 
and  vertical  velocity  were  approximated  from  precalculated  values 
instead  of  observed  values  of  instruments,  the  RMS  errors  were  not 
estimated. 
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6.  RESULTS 

a.  Statl sties 

Samples  of  the  charts  depicting  prestorm  fields  are  contained  In 
Appendix  A.  Figures  A-la  through  A-lk  Illustrate  the  fields  associated 
with  the  15  August  1982  MCC.  Figures  A-2a  through  A-2k  are  similar  but 
for  the  9  August  1978  non-MCC  storm.  Table  8  summarizes  the  Informa¬ 
tion  obtained  from  the  eighteen  storm  cases.  It  lists  the  range  and 
mean  value  of  each  variable  for  the  two  storm  categories. 

By  comparing  these  MCC  and  non-MCC  values,  several  differences 
can  be  seen.  Looking  at  the  range  column  first,  note  that  considerable 
overlap  occurs  with  each  variable.  While  the  most  favorable  values  for 
strong  convection  are  associated  predominantly  with  the  MCC,  this  Is 
not  always  the  case.  The  non-MCC  category  has  a  higher  value  of  850  mb 
water  vapor  flux  convergence  and  a  higher  Totals  Index  value.  Its 
maximum  value  for  surface-to-500  mb  average  relative  humidity  equals 
that  of  the  MCC  column.  On  the  opposite  end  of  the  scale,  those 
values  least  favorable  for  convection  are  generally  associated  with  the 
non-MCCs,  vortlclty  advectlon  and  forecast  vertical  velocity  being  the 
only  exceptions. 

Referring  to  the  mean  value  column,  a  few  more  relationships  are 
evident.  While  there  are  no  sign  differences  between  the  correspond¬ 
ing  MCC  and  non-MCC  values,  the  MCC  mean  value  In  every  category  is 
more  conducive  to  convection.  For  the  vapor  variables,  two  important 
differences  should  be  noted.  First,  the  water  vapor  content  of  the  MCC 
prestorm  environment  Is  23  percent  greater  than  that  for  the  non-MCC 


Table  8.  Variable  range  and  mean  value  for  MCC  and  non-MCC  cases. 
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number  of  height  contour  and  vorticity 
isoline  intersections) 
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environment.  This  difference  is  reflected  by  both  the  850  mb  mixing 
ratio  and  preclpitable  water.  Second,  the  MCC  low-level  water  vapor 
processes  appear  much  more  vigorous.  Vapor  advection  is  more  than 
three  times  greater  while  vapor  flux  convergence  is  almost  twice  as 
great.  The  convergence  value  of  .46  g  kg“*h“*  compares  fairly  well 
with  Maddox's  (1981)  prestorm  value  of  0.32  g  kg'*  h~*.  A  signifi¬ 
cant  difference  exists  between  the  advection  values,  however.  Maddox 
obtained  a  value  of  about  -0.03  g  kg“*  h”*  compared  to  the  +0.20  g 
kg"*  h"1  obtained  in  this  study.  A  significant  part  of  this 
difference  may  be  due  to  either  sample  size  or  differing  evaluation 
techniques.  Maddox  determined  the  average  value  over  an  area  larger 
than  the  storms  IR  satellite  image.  Maximum  values  occurring  within 
the  smaller  radar  echo  region  were  used  for  this  study.  The  final 
observation  concerns  thickness  advection.  Mean  low-level  warm 
advection  is  more  than  three  times  greater  for  the  MCC  cases.  While 
upper-level  cold  advection  is  greater  also,  the  difference  is  much 
smaller.  Compared  to  earlier  studies,  the  low-level  warm  advection 
value  of  +2.5  gpm  h“*  is  slightly  low.  Maddox  (1983)  reported  an 
approximate  value  of  +3.5  gpm  h“*  in  his  investigation  of  ten  MCCs. 
Hales  (1982),  in  an  investigation  of  seventeen  tornado  producing 
storms,  reported  an  average  value  of  +3.0  gpm  h-*. 

b.  Statistical  Comparison  of  Means 

Results  from  the  statistical  comparison  of  MCC  to  non-MCC  mean 
values  are  shown  in  Table  9.  For  each  variable,  the  a  level  is  listed 
at  which  the  two  mean  values  are  determined  to  be  significantly 


Table  9.  Alpha  (a)  levels  at  which  the  corresponding  MCC  and  non-MCC 
mean  values  were  determined  to  be  significantly  different.  Only 
levels  between  0  and  0.5  were  considered. 
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Parameter  Level 


850  mb  mixing  ratio  .10 
850  mb  advection  of  water  vapor  .05 
850  mb  flux  divergence  of  water  vapor  .40 
Surface-to-500  mb  average  relative  humidity  None 
Preci pi  table  water  .10 
Totals  Index  .20 
Total  Energy  Index  .10 
1000-700  mb  thickness  advection  .05 
500-300  mb  thickness  advection  None 
500  mb  vorticity  advection  None 
700  mb  vertical  velocity  forecast  .20 


different.  These  a  levels  express  the  probability  of  erroneously  con¬ 
cluding  that  the  two  values  come  from  different  populations.  There¬ 
fore,  those  variables  having  large  a  levels  should  be  of  little  value 
in  distinguishing  between  MCC  and  non-MCC  prestorm  environments. 
Statistical  results  imply  that  surface-to-500  mb  average  relative 
humidity,  500-300  mb  thickness  advection,  and  vorticity  advection  are 
of  no  value.  When  initially  comparing  mean  values,  vapor  flux  conver¬ 
gence  appeared  promising  as  a  distinguishing  variable.  The  test  for 
significant  difference  does  not  support  this  contention,  however.  The 
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large  variability  of  vapor  flux  convergence  makes  it  less  effective. 
Vertical  motion  forecasts  and  the  Totals  Index  show  some  value  with  an 
a  level  of  .20.  Of  the  two  stability  indices,  however,  the  Total 
Energy  Index  appears  more  useful  with  an  a  level  of  .10.  Other 
variables  having  an  a  level  of  .10  include  the  850  mb  mixing  ratio  and 
precipi table  water.  The  two  variables  demonstrating  the  greatest 
potential  for  distinguishing  between  MCC  and  non-MCC  environments  are 
low-level  thickness  and  low-level  water  vapor  advection.  Both  were 
determined  significant  at  the  .05  level. 

c.  Forecast  Applications 

The  application  of  these  results  to  the  realm  of  forecasting 
should  be  viewed  with  caution  due  to  the  small  sample  size.  Regard¬ 
less,  possible  applications  should  be  identified.  In  an  attempt  to 
distinguish  between  the  MCC  and  non-MCC  prestorm  environments,  scatter 
diagrams  were  used.  Figure  7  illustrates  one  of  these  diagrams.  For 
this  diagram,  low-level  thickness  advection  is  plotted  against  upper- 
level  thickness  advection  for  all  eighteen  storms.  The  nine  MCCs  are 
numbered  and  represented  with  circles.  The  nine  non-MCCs  are  numbered 
also  but  represented  with  triangles.  The  storm  numbers  correspond  to 
the  numbers  listed  in  Table  5.  While  the  non-MCCs  tend  to  cluster, 
the  intermingling  of  several  MCCs  prevents  a  clear  separation.  This 
poor  separation  is  not  surprising,  however,  since  upper-level  thickness 
advection  was  determined  to  be  an  ineffective  indicator  for  distin¬ 
guishing  between  the  two  storm  types. 

Since  low-level  thickness  advection  and  water  vapor  advection 


500-300  MB 
THICKNESS 
ADVECTION 
BY  400  MB 


1000-700  MB  THICKNESS  ADVECTION 
BY  850  MB  WIND 
(gpm  h-D 

Fig.  7.  Scatter  diagram  plot  of  1000-700  mb  and  500-300  mb 
thickness  advection  values  for  MCC  (circles)  and  non-MCC  (triangles) 
storms.  Storm  numbers  correspond  to  numbers  in  Table  5. 
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had  the  lowest  a  levels,  a  plot  of  these  two  variables  should  provide 
the  best  separation.  This  scatter  diagram  is  depicted  in  Fig.  8  and 
does  provide  the  best  separation  when  compared  to  other  possible 
combinations.  With  the  exception  of  MCC  numbers  3  and  5,  a  separation 
is  evident. 

In  hopes  of  enhancing  this  separation,  individual  variables  were 
combined  in  a  physically  realistic  manner  and  replotted.  This  approach 
is  best  illustrated  by  returning  to  Fig.  7.  Excluding  MCC  number  3,  a 
relation  between  MCC  low-level  and  upper-level  advection  was  noted. 
Those  MCCs  characterized  by  strong  low-level  warm  advection  showed 
little  advection  aloft.  However,  those  MCCs  with  weak  low-level  warm 
advection  had  significant  compensating  values  of  cold  advection  aloft. 
The  non-MCCs  did  not  demonstrate  a  similar  relationship.  Thus,  a  term 
including  both  low  and  upper-level  advection  might  provide  a  better 
separation  of  storm  types.  This  idea  was  tested  in  the  following 
manner.  A  single  term  was  formed  by  subtracting  the  500-300  mb  thick¬ 
ness  advection  value  from  the  corresponding  1000-700  mb  thickness 
advection  value.  Although  equivalent  rates  of  thickness  advection 
within  the  two  layers  produce  different  rates  of  mean  virtual  tempera¬ 
ture  advection,  the  difference  was  not  considered  large  enough  to  make 
weighting  necessary.  The  resulting  term,  named  the  Differential 
Advection  Index  (DAI),  provides  a  measure  of  the  change  in  stability 
with  time.  Positive  values  indicate  that  the  lapse  rate  is  destabiliz¬ 
ing  with  time;  the  larger  the  value,  the  greater  the  rate  of  change 
will  be.  After  determining  the  DAI  value  for  each  storm,  means  were 
computed  for  the  MCCs  and  non-MCCs.  Applying  the  Student  t  test,  the 
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Fig.  8.  Scatter  diagram  plot  of  1000-700  mb  thickness  advection  and  850  mb  water  vapor 
advection  for  MCC  (circles)  and  non-MCC  (triangles)  storms.  Storm  numbers  correspond  to  numbers 
in  Table  5.  When  used  as  a  forecast  nomogram,  the  area  above  the  dashed  line  represents 
conditions  favorable  for  MCC  formation. 
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means  were  determined  to  be  significantly  different  at  the  .10  a  level. 
Thus,  the  combination  process  greatly  Increased  the  value  of  upper- 
level  advectlon  as  a  distinguishing  factor.  Figure  9  shows  the  storm 
separation  achieved  by  plotting  the  850  mb  mixing  ratio  against  the 
DAI.  Compared  to  Fig.  8,  this  scatter  diagram  does  a  better  job  of 
separating  MCC  numbers  3  and  5  from  the  non-MCC  storms.  Non-MCC 
number  4,  however,  is  now  located  in  the  midst  of  the  MCC  region.  It 
is  worth  noting  that  number  4  was  a  large  thunderstorm  system  that 
developed  in  the  same  region  in  which  an  MCC  had  occurred  the  previous 
day.  The  storm's  environment  possibly  was  still  quite  similar  to  that 
of  the  MCC. 

Since  low-level  thickness  and  water  vapor  advection  had  provided 
fairly  good  separation  when  plotted  against  one  another,  an  effort  was 
made  to  combine  these  two  variables.  A  clue  to  combining  these  two 
terms  was  provided  by  Darkow's  development  of  the  Total  Energy  Index. 
Mis  index  is  based  upon  the  difference  in  static  energy  at  the  500  mb 
and  850  mb  levels.  The  energy  totals  at  each  level  are  determined  by 
calculating  the  contributions  of  sensible  heat,  latent  heat,  and 
potential  energy.  If  potential  energy  is  ignored,  a  method  is  thus 
available  for  combining  the  thickness  and  water  vapor  advection 
variables.  Using  advection,  however,  gives  an  estimate  of  the  rate  of 
change  of  energy  rather  than  an  actual  energy  total.  This  method  is 
described  in  Appendix  B.  Thus,  the  energy  rate  of  change  at  850  mb 
was  computed  for  each  storm,  averaged  for  the  two  storm  types,  and 
then  tested  for  its  ability  to  distinguish  between  MCC  and  non-MCC 
storm  environments.  The  means  were  determined  to  be  significantly 
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Fig.  9.  Scatter  diagram  plot  of  850  mb  mixing  ratio  and  Differential  Advection  Index  for 
MCC  (circles)  and  non-MCC  (triangles)  storms.  Storm  numbers  correspond  to  numbers  In  Table  5. 
When  used  as  a  forecast  nomogram,  the  area  above  the  dashed  line  represents  conditions  favorable 
for  MCC  formations. 
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different  at  the  .01  a  level. 

Additional  scatter  diagrams  were  drawn  utilizing  the  850  mb 
energy  rate  of  change.  Figure  10  shows  the  850  mb  mixing  ratio  plotted 
against  the  energy  rate  of  change.  While  grouping  is  good,  differences 
between  MCC  number  5  and  non-MCC  number  4  remain  unresolved.  Figure 
11  shows  a  plot  of  the  energy  rate  of  change  against  the  adjusted 
Total  Energy  Index.  Using  the  energy  rate  of  change  values  previously 
calculated  for  each  storm,  850  mb  static  energy  totals  were  computed 
for  the  time  0000  GMT  +7  h.  In  turn,  these  850  mb  energy  totals  were 
used  to  compute  an  adjusted  Total  Energy  Index  value  for  each  storm. 

The  projected  stability  values  were  valid  closer  to  the  MCCs'  formation 
times  and,  thus,  hopefully  more  representative  of  the  formation  envi¬ 
ronment.  Regardless,  scatter  diagram  results  remained  essentially  the 
same.  Figure  12  depicts  the  final  scatter  diagram.  The  energy  rate 
of  change  has  been  plotted  against  the  DAI.  Grouping  is  quite  good 
but,  again,  complete  separation  is  not  attained.  Table  10  provides  a 
summary  of  the  new  variables  derived.  In  addition,  it  depicts  the 
corresponding  MCC  and  non-MCC  mean  values  and  the  a  levels  at  which 
these  mean  values  are  determined  to  be  significantly  different. 

Based  upon  the  scatter  diagram  results,  Figs.  8,  9,  10,  and  12 
are  recommended  for  use  as  experimental  MCC  forecast  nomograms.  For 
each  one,  the  region  above  the  dashed  line  represents  environmental 
conditions  favorable  for  MCC  development.  No  single  nomogram  incor¬ 
porates  all  of  the  variables  found  useful  in  distinguishing  between 

MCC  and  non-MCC  prestorm  environments.  Each  provides  a  different 

i 

perspective.  Therefore,  probably  all  of  the  graphs  should  be  used 
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10.  Scatter  diagram  plot  of  850  mb  mixing  ratio  and  850  mb  energy  change  rate  for  MCC 
and  non-MCC  (triangles)  storms.  Storm  numbers  correspond  to  numbers  in  Table  5.  When 
forecast  nomogram,  the  area  above  the  dashed  line  represents  conditions  favorable  for 


11.  Scatter  diagram  plot  of  850  mb  energy  change  rate  and  adjusted  Total  Energy  Index 
rcles)  and  non-MCC  (triangles)  storms.  Storm  numbers  correspond  to  numbers  in  Table  5 
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Fig.  12.  Scatter  diagram  plot  of  850  mb  energy  change  rate  and  Differential  Advection  Index 
for  MCC  (circles)  and  non-MCC  (triangles)  storms.  Storm  numbers  correspond  to  numbers  in  Table 
5.  When  used  as  a  forecast  nomogram,  the  area  above  the  dashed  line  represents  conditions 
favorable  for  MCC  formation. 
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Table  10.  The  new  variables  formed,  their  mean  MCC  and  non-MCC 
values,  and  the  alpha  (a)  level  at  which  each  pair  was  determined  to 
be  statistically  different. 


PARAMETER 

MEAN 

MCC 

VALUE 

NON-MCC 

a  LEVEL 

Differential  Advection  Index 
(gpm  h-D 

4.6 

2.1 

.10 

850  mb  Energy  Rate  of  Change 
(cal  g-l  h~l) 

.18 

.05 

.01 

Adjusted  Total  Energy  Index 
(cal  g_1) 

-5.0 

-2.9 

.01 

when  evaluating  MCC  formation  potential;  this  procedure  is  especially 
important  for  borderline  storm  environments.  These  graphs  should  be 
used  for  only  short  term  forecasts  since  they  assume  adiabatic 
conditions  and  perfect  advection. 
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7.  CONCLUSIONS 

Statistical  results  suggest  that  certain  atmospheric  meteorologi¬ 
cal  variables  can  be  used  to  distinguish  between  MCC  and  non-MCC 
prestorm  environments.  Variables  of  statistical  value  include  the  850 
mb  mixing  ratio,  850  mb  water  vapor  advection,  precipitable  water,  the 
Total  Energy  Index,  and  1000-700  mb  thickness  advection.  The  MCC  and 
non-MCC  mean  values  for  each  of  these  variables  were  determined  to  be 
significantly  different  at  n  a  level  of  .10  or  less. 

The  scatter  diagram  plots  reinforce  the  statistical  findings. 

Of  the  eleven  original  variables,  low-level  thickness  advection  and 
low-level  water  vapor  advection  provide  the  best  graphical  separation 
of  MCC  and  non-MCC  storms.  New  indicators  developed  by  combining 
variables  also  show  skill  in  graphically  separating  the  two  storm 
types.  While  upper-level  advection  by  itself  is  a  poor  indicator,  its 
combination  with  low-level  advection  provides  a  more  complete  measure 
of  differential  advection  and,  thus,  is  more  useful.  Low-level 
thickness  and  water  vapor  advection  have  been  combined  to  provide  an 
estimate  of  the  energy  change  rate  at  850  mb.  Of  all  the  variables 
investigated,  this  energy  rate  of  change  term  provides  the  best 
distinction  between  MCC  and  non-MCC  prestorm  environments.  This 
result  implies  that  abnormally  large  energy  increases  within  the 
low-levels  are  required  for  MCC  development. 

MCC  and  non-MCC  storm  separation  is  apparent  on  the  scatter 
diagrams.  While  complete  separation  is  never  attained,  the  results 
are  encouraging,  especially  considering  the  difficulty  of  the  task. 

One  cannot  expect  synoptic  data  to  provide  a  perfect  distinction 


between  the  environments  of  two  mesoscale  events  so  similar  In  nature. 

While  Investigative  results  are  encouraging,  they  cannot  be 
accepted  with  a  large  degree  of  confidence.  In  general,  differences 
between  the  MCC  and  non-MCC  mean  values  are  only  slightly  larger  than 
the  associated  RMS  errors.  Furthermore,  the  small  number  of  cases 
examined  reduces  the  reliability  of  the  statistical  test  applied.  A 
final  consideration  concerns  storm  location.  Seven  of  the  nine  MCCs 
studied  formed  In  either  eastern  Kansas  or  Missouri.  Thus,  mean  values 
obtained  for  the  MCCs  may  have  limited  application  elsewhere, 
especially  for  the  drier  western  plains  region. 

Additional  Investigations  considering  a  larger  number  of  storms 
are  needed  to  verify  the  results  obtained  here.  A  greater  number  of 
case  studies  would  allow  a  more  reliable  statistical  analysis.  If 
results  continue  to  be  encouraging,  a  sound  basis  will  exist  for  the 
development  of  operational  forecast  procedures  and  techniques.  While 
synoptic  data  do  not  provide  sufficient  Information  to  pinpoint  the 
location  and  timing  of  severe  weather  events,  they  have  proven  useful 
In  Identifying  those  areas  most  susceptible  to  such  weather.  Use  of 
synoptic  data  to  identify  those  environments  favorable  for  MCC  develop¬ 
ment  appears  feasible.  Such  guidance  could  have  a  significant  impact 
upon  regional  forecast  preparation,  especially  when  considering  the 
potential  for  heavy  rains,  flooding,  hail  and  tornadoes. 
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APPENDIX  A 

CHARTS  FOR  THE  15  AUGUST  1982  MCC 
AND  9  AUGUST  1978  NON-MCC 

Grid  values  of  the  first  nine  variables  listed  in  Table  2  were 
computed  for  each  MCC  and  non-MCC  case;  isopleths  were  then  drawn  by 
hand.  NMC  500  mb  Height/Vortici ty  Analysis  and  LFM  vertical  velocity 
forecast  panels  were  used  to  represent  the  two  final  variables.  Samp¬ 
les  of  these  charts  and  NMC  preducts  are  contained  in  this  appendix. 
Figures  A-la  through  A-lk  illustrate  the  fields  associated  with  the  15 
August  1982  MCC.  Figures  A-2a  through  A-2k  illustrate  the  fields 
associated  with  the  9  August  1978  non-MCC. 
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Fig.  A-la.  850  mb  mixing  ratio  (g  kg_M  for  0000  GMT  15 
August  1982.  Stippled  region  depicts  the  MCC  radar  echo  at  0835  GMT. 
Radar  echoes  not  a  part  of  the  MCC  are  not  shown. 
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Fig.  A-lb.  850  mb  water  vapor  advection  (xl0~2  g  kg-*  h~l) 
for  0000  GMT  15  August  1982.  Solid  lines  are  positive  and  dashed 
lines  are  negative.  Stippled  region  depicts  the  MCC  radar  echo  at 
0835  GMT. 
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Fig.  A-lc.  850  mb  flux  divergence  of  water  vapor  (xl0~2  g 
kg- 1  h“l)  for  0000  GMT  15  August  1982.  Solid  lines  are  positive 
and  dashed  lines  are  negative.  Stippled  region  depicts  the  MCC  radar 
echo  at  0835  GMT. 


Fig.  A-ld.  Surface-to-500  mb  average  relative  humidity  (%)  for 
0000  GMT  15  August  1982.  Stippled  region  depicts  the  MCC  radar  echo 
at  0835  GMT. 


Fig.  A-lg.  Total  Energy  Index  (cal  g~l)  for  0000  GMT  15  August 
1982.  Stippled  region  depicts  the  MCC  radar  echo  at  0835  GMT. 
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000-700  mb  thickness  advection  (gpm  Ir*)  by  850 
5  August  1982.  Solid  lines  are  positive  and  das 
Stippled  region  depicts  the  MCC  radar  echo  at 
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Fig.  A-li.  500-300  mb  thickness  advection  (gpm  h-*)  by  400  mb 
wind  for  0000  GMT  15  August  1982.  Solid  lines  are  positive  and  dashed 
lines  are  negative.  Stippled  region  depicts  the  MCC  radar  echo  at 
0835  GMT. 
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Fig.  A-lj.  NMC  500  mb  Height/Vorticity  Analysis  Chart  for  0000 
GMT  15  August  1982.  Units  are  gpm  for  height  (solid  lines)  and 
(xlO'5  s-1)  for  vorticity  (dashed  lines).  The  height/vorticity 
field  within  the  box  was  used  to  estimate  the  type  and  intensity  of 
vorticity  advection  into  the  MCC  region. 


vertical  velocity  are  dashed. 


Fig.  £-2h.  1000-700  mb  thickness  advection  (gpm  h-*)  by  850  mb 
wind  for  0000  GMT  9  August  1978.  Solid  lines  are  positive  and  dashed 
lines  are  negative.  Stippled  region  depicts  the  non-MFC  radar  echo  at 
0435  GMT. 
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Iso-lines  of  vertical  velocity  are  dashed. 
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APPENDIX  B 

METHOD  USED  TO  ESTIMATE  THE  850  MB  ENERGY 
RATE  OF  CHANGE  DUE  TO  THICKNESS  AND  WATER  VAPOR  ADVECTION 

The  static  energy  (Es)  of  a  unit  mass  of  air  is  expressed  as 

Es  -  CpT  t  Lq  *  gZ.  (B1) 

It  is  the  sum  of  specific  enthalpy  (cpT),  latent  heat  (Lq),  and 
potential  energy  (gZ).  Darkow  (1968)  approximated  this  quantity 
using  (B2) 

Es  “  cpT  +  Low  +  gZ  (B2) 

where  LQ  is  a  constant  latent  heat  of  condensation  and  w  is  the  mixing 
ratio.  The  energy  rate  of  change  is  expressed  by 


dT 

It 


+  L, 


dw 

It 


+  g 


dZ 

IT 


(B3) 


For  this  investigation,  dEs  is  defined  as  the  energy  rate  of  change 

dt 

due  solely  to  changes  In  sensible  and  latent  heat.  Thus, 


dE 

It' 


s  = 


CP  IT  + 


dw 

IF 


(B4) 


Furthermore,  If  only  advective  changes  are  considered,  (B4)  becomes 


Equation  (B5)  provides  a  means  for  combining  low-level  thickness  and 
water  vapor  advection  into  one  physically  realistic  term.  This  proce¬ 
dure  was  applied  to  each  MCC  and  non-MCC  case  to  obtain  estimates  of 
the  850  mb  energy  rate  of  change  due  to  advection. 

Before  (B5)  could  be  used,  however,  thickness  advection  values 
had  to  be  converted  to  advection  values  of  mean  virtual  temperature 

( T* )  using  the  relationship 

(at*\  |  adv  ( B6 } 

Ut  jadv  '  K  ln  £l 
Pz 

« 

where  K  is  the  hypsometric  constant.  The  expression  AV(A4>)  represents 
the  vertical  change  in  thickness  between  pressure  levels  p!  and  p2. 
Simplified  for  calculations,  (B5)  becomes 

■  •*«  (if)  ad,  *  -586  (u)  adv  <B7) 

\ At  /  adv 

where  units  for  /Al  \  an<^  (—  1  are  °C  h_l  and  g  kg~^  h"l, 

\ a t  /adv  V At /  adv 

respectively.  L0  was  selected  based  upon  a  typical  850  mb  temperature 
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of  20°C.  Units  for  /aE£\  •  the  energy  rate  of  change,  are 

\ At  /  adv 

cal  g_1  h_1. 
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